Coexistence of s-wave Superconductivity and Antiferromagnetism 
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We study the phase diagram of a new model that exhibits a first order transition between s-wave 
superconducting and antiferromagnetic phases. The model, a generalized Hubbard model augmented 
with competing spin-spin and pair-pair interactions, was investigated using the projector Quantum 
Monte Carlo method. Upon varying the Hubbard U from attractive to repulsive we find a first order 
phase transition between superconducting and antiferromagnetic states. 
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Experimental phase diagrams with superconducting 
and magnetic phases in close proximity continue to fas- 
cinate since the broken symmetry states are incompat- 
ible. In particular, a canonical two-dimensional ex- 
ample is the layered organic superconductor k-(BEDT- 
TTF)2Cu[N(CN)2]Cl which shows a first order transi- 
tion between antiferromagnetic (AF) insulating and su- 
perconducting (SC) phases upon varying the pressure. |l] 
Debate continues about the nature of the superconduct- 
ing state. Kanoda et al. argue that NMR data are con- 
sistent with an unconventional superconducting gap with 
nodal lines. Q Conversely, many experimental studies are 
rather consistent with s-wave, phonon mediated super- 
conductivity. In particular, Kini et al. could demonstrate 
a BCS-like mass isotope effect in k-S.|^ The same mate- 
rial exhibits a strong superconductivity induced (acous- 
tic) phonon renormalization. Q Nodes in the pair or- 
der parameter are also incompatible with specific heat 
measurements. Q 

In order to model the scenario of a phase transition 
between an AF and an s-wave SC on a half-filled two di- 
mensional lattice we add to the standard Hubbard model 
a pair hopping term: 



where i labels the sites of a square lattice and the first 
sum runs over nearest neighbor bonds. The fermionic 
spin is given by = 1 ^ c\ ^(Ts,tC^^ and a denotes the 
Pauli matrices. For a dimer the pair hopping term arises 
exactly when we assume a Su-Schrieffer-Heeger electron- 
phonon interaction |^ and integrate out the phonon 
in the antiadiabatic limit. We introduce a particle- 
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pairing operators as 0] r/i" 



VS:fP , the pair creation op- 
erator 77-1 (—l)*c-!. cl and 7?5 = (ri^-f-f 77^ I — 1)/2 is the 
total charge. The pair hopping term may be recast 
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simulations, we have used the projector auxiliary field 
QMC algorithm which is well suited for the study of 
ground state properties Simulations do not suffer 
from a sign problem in the region of attractive U < 
with the additional constraint U < — 8|y|. This is the 
shaded region in Fig. Using particle-hole symmetry 
the region of repulsive J7 at = is also accessible. 
Two SU{2)s,ri symmetry groups are generated by the 



spin algebra 
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2(5jjS'| and its particle- hole 



equivalent for the rj operators. Note that the commu- 
tator 5^, ff:^ = because S"?]^ = 0. The hopping terms 
Ht and Ht^, are known to be invariant under spin rotation 
and also commute with the discrete symmetry V. This 
in turn implies that 77^ = VSfP also commutes with both 
hopping terms. The symmetries of the Hubbard term 
become apparent after putting it into the form 
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hole transformation V with properties Pcj 



We can read off the SU{2)s ® SU{2)n symmetry and ob- 
serve that the Hubbard term changes sign under V trans- 
formation. 10] In Fig. Qthe various symmetries for differ- 
ent parameter regions are displayed. Bars denote symme- 
tries which we numerically find to be spontaneously bro- 
ken in the thermodynamic limit. The point U = V = Q 
has the full symmetry: SU{2)s® SU (2)^ (g) ^2,p_f/ where 
Z2.PH denotes the discrete particle-hole symmetry. Here 
two scenarios are possible: i) Long range order for both 
AF and SC correlations since one implies the other 
through Z2,pH symmetry or ii) The competition between 
the two broken symmetry states leads to a disordered 
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FIG. 1: The phase diagram for the full Hamiltonian. The 
shaded region indicates the parameters accessible to sign-free 
simulations. For every region in the phase diagram the full 
symmetry of the Hamiltonian is indicated. Broken symme- 
tries are marked by bars. 



state. We find that the system realizes the first possi- 
bility of coexistence of s-wave SC and AF. Switching on 
U breaks the particle-hole symmetry whereas V reduces 
the SU{2)^ U{1) (g) Z2,r,. 

In the following, we will first discuss the phase transi- 
tion at constant negative J7 as a function of V and then 
turn our attention to the phase transition at V = and 
varying U. 

In the U ^ ~oo limit the model is equivalent to hard- 
core bosons with repulsive nearest neighbor interactions. 
At half-filling this system experiences a transition from 
SC to charge-density wave states (CDW) going across 
V — 0. 1^ The transition point is marked by the higher 
symmetry SU (2)^ and broken symmetry states may be 
labeled according to the direction of the magnetization. 
Therefore it is possible to find broken symmetry states 
which show a coexistence of CDW and SC order param- 
eters. One can view this as a consequence of merging the 
Ising and the U{1) symmetries into a continuous group. 

The same behavior is still present for finite attractive 
U. In particular our simulations show (see Figs, and 
Pd) that single particle excitations are always gapped so 
we can argue that the system still renormalizes to the 
bosonic model. As we will see, crossing the ^ = line at 
T = 0, we observe a divergence of the correlation length 
(Fig. 1^) which points to a second order phase transition. 
Furthermore, we find no sign of a discontinuity in the first 
derivative of the free energy {l/N)dF/dV = ^2(n^j^-) 
(see Fig. ^) . At the same time we observe a jump for the 
CDW and SC order parameters at F = 0. To understand 
this transition, we split it up into three parts: a) V ^ 0~, 
b) V = 0" ^ 0+ and c) V ^ . Part b) is responsible 
for the jump in the order parameters. Let us introduce 



^-ising = hm lim *(F) 

y^O- N^oo 

^fxY = hm lim '^(V). 
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The ground state in the Ising regime ^'ising has CDW or- 
der and the SC order parameter is zero. Jumping to posi- 
tive V and ^xY the situation reverses. Thus the observed 
jump in the order parameters is the simple consequence 
of forcing the transition from one broken symmetry state 
to the other: Rising ^xy- Next we look at the critical 
behavior as we approach the phase transition from either 
side. Linear spin wave theory (LSWT) predicts for both 
transitions a) and c) the same square root divergence for 
the respective transverse correlation function. In case c) 
the equal time transverse correlation function is {n^n_^) 

where = {\/^/N)^pe^'^^np, np = + np^^ and 

Q — (tt, tt). According to LSWT this density-density cor- 
relation diverges like {n^n_^) — %S/ ^Jvjtp as F ^ 0+. 
The real part of the uniform susceptibility 



X' „(^ = Q'^ = 0) = 



dr 



5)ir) (5) 



picks up the inverse of the gap A — 8y^2Vtp and in 
LSWT diverges hke Xn,niQ' 0) = ^/T/2S/V. Asymptoti- 
cally the correlation function {npiQ^ (r) decays exponen- 
tially in space and time thus the integrated correlation 
functions are proportional to the correlation length and 
gap 
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Data for both Xn.niQ^ 0) and {n^n^ ^)(t = 0) are shown 
in Fig. ^ where the dotted line plots the LSWT result 
multiplied with a single constant to account for the re- 
duced moment. 

The ordered state in the XY regime has one Gold- 
stone mode and the transverse correlations (?t-q't-_q) are 
gapped. Approaching the critical point the gap closes 
and the {n(^n_Q) correlation length diverges. Thus crit- 
ical behavior reflects the softening of the observed mode 
which turns into a Goldstone mode at the point with 
higher symmetry. Ultimately it is the softening of the 
transverse mode which is responsible for the observed 
critical behavior. 

Following the V — line as [/ ^ the bosonic LRO 
decreases only slightly leading to a finite limit at C/ = 0+ 
which is shown in Fig. |^. Crossing the U = point we 
encounter again a phase transition which is linked to a 
point with higher symmetry: the discrete particle-hole 
symmetry. Using the Z2.PH symmetry, one can map this 
bosonic LRO to an equivalent spin LRO for the repulsive 
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FIG. 2; a) Uniform susceptibility Xn,n(O,0) oc ^^A"^ b) 
oc The dotted line is proportional to results 



from LSWT. 
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FIG. 3: a) Density-density correlations 2 {nin:^_^_g) — 
— l/N{dF/dV) where = rii^j + The derivative of the 
free energy is continuous, b) The quasi-particle gap A^p was 
obtained by fitting the tail of the imaginary time correlation 
function G(r). 



U case at y = 0. Precisely at J7 = particle hole sym- 
metry implies that bose correlation functions are equal 
to their spin counterparts 



(8) 



This explains a jump by a factor two for the density- 
density correlations as one moves from the symmetric 
point to a finite U which is visible in Fig. ^. The emer- 
gence of a Goldstone mode above was associated with the 
critical behavior of the phase transition. As we are now 
approaching the U — point the only continuous sym- 
metry available which is not associated with the bosonic 
LRO is the spin symmetry. In Fig. ^ we compute the cor- 
relation length for {S^S^' -). This quantity shows no sign 
of divergence at the phase transition. Contrary to what 
happened in the SC to CDW transition we may classify 



FIG. 4: The density- density correlation function 
1/N {n{L/2)n{0)) at maximum separation (L/2,L/2). 
The correlation function jumps by a factor two once a finite 
Hubbard U term breaks the Z2,ph symmetry. 
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FIG. 5: Values for the uniform susceptibility x's^.s^i^ ~ 
Q,uj — 0) cx ^^^^ (lower curve) are well below the correla- 



tion function {51 SI-) (A; = Q) oc C 



this transition as a first order level crossing transition. 
From Fig. || we equally see that the derivative of the 
free energy with respect to U is discontinuous. The free 
energy as a function of double occupancy has a plateau 
which marks a region of phase coexistence. 

For the point U = V — Owe find coexistence of the spin 
and pair order parameters. After finite size extrapolation 
the results of Fig. ^ are consistent with a finite value of 
the magnetization rrig up to weak coupling (tp/t << 1). 
At the same time the superconducting order parameter 
A is also nonzero and the measured value A = y/2ms 
agrees with Eq. |[ The observed coexistence of AF and 
SC order is due to the special particle-hole symmetry 
Z2,PH at the point U = which renders the energies of 
both the AF and SC states degenerate. 

In order to understand how the SU (2)5 (g) SU (2)^ sym- 
metries are spontaneously broken we also need to discuss 
whether the discrete particle-hole symmetry is broken 
and how this affects the spontaneous breaking of the con- 
tinuous symmetries. We introduce the following pseudo- 
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FIG. 6: a) The "Hubbard" correlation (n|nx) = 

l/N{dF/dU) has a jump at the transition point U — 0. The 
first order transition is signaled by a discontinuity of the first 
derivative of the free energy, b) The quasi-particle gap Aqp. 
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FIG. 8: The l/iV(5g5g)(fc = 0) correlations show that at 
least up to T < l.St particle-hole symmetry is broken. 
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FIG. 7: The staggered moment rris which coexists with the 
superconducting order parameter A = msV^ plotted as a 
function of tp {U/t = V/t = 0). The inset shows finite size 



extrapolation for tp — i and tp 



0.5. 



In conclusion we have considered a model with com- 
peting interactions. The main ingredients in the model 
are on the one hand the Heisenberg and Hubbard terms 
which lead to a localization of single electrons thereby 
producing antiferromagnetic insulating states and on the 
other hand, pair hopping and nearest neighbor attrac- 
tive terms which lead to on-site s-wave superconducting 
states. Based on the symmetries of the model as well as 
quantum Monte Carlo simulations we have shown that 
the model has a variety of phase transitions. In particu- 
lar we have found a first order phase transition between 
an s-wave superconducting state and antiferromagnetic 
Mott insulator at a high symmetry point {U = V = 0) 
where the two phases coexist. 
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{S) = {VSV) = - {S) 
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which has to be zero as long as the Z2.PH symmetry 
is not broken as assumed in Eq. (|lO|). A positive 
(negative) value for {S) indicates a majority of spins 
(pairs) in the system. Plotting the uniform correlations 
1/N (^Sj:Sj:) (k = 0) as a function of temperature, one 
sees in Fig. ^ that below a temperature Tc ~ \M particle- 
hole symmetry is indeed broken. Hence below Tc and in 
the thermodynamic limit the Z2^ph particle-hole sym- 
metry is broken and depending on the orientation of the 
pseudo-spin, 5^, the SC or AF ground state will be cho- 
sen. 
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